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I. INTRODUCTION 



The search for new physics beyond the Standard Model (SM), especially supersymme- 
try (SUSY), is one of the objectives at the CERN Large Hadron Collider (LHC). Many 
calculations have been carried out based on the Minimal Supersymmetric Standard Model 
(MSSM), a version of SUSY. Phenomenologically SUSY predicts many new particles; e.g., 
the superpartners of the SM particles. Specifically, in the MSSM, there are squarks, gluino, 
sleptons, charginos, neutralinos and more Higgs bosons in addition to the SM particles. Be- 
sides squarks and gluinos, perhaps the most interesting new particles are the four neutralinos 
Xiii = 1,2,3,4) and the two charginos xfU = 1,2), which are the mass eigenstates of the 
superpartners of the Higgs and gauge bosons, since the lightest chargino xf and the two 
lightest neutralinos {xi, X2) can be lighter than the squarks and the gluino in most of the 
parameter space. In most of the MSSM parameter regions the associated production of xfx2 
is the main source of trilepton events. In Refs.[l|, Q] the trilepton signal was investigated 
for the Fermilab Tevatron in the Minimal Supergravity Model (mSUGRA) at leading order 
(LO). If the leptonic decays of xfxt are the dominant decay modes the signal to background 
ratio can be quite large after suitable cuts. Because the trileptonsignal is also quite sensitive 
to the SUSY parameters, it is potentially also a sensitive probe of the SUSY parameters. In 
fact, the trilepton signal from xfx2 ^'^^ being searched for by the DO Collaboration 3| at 
the Tevatron. So far no excess has been observed above the expected SM background, but 
the results have been used to constrain the masses. A plan for searching for the trilepton 
signal from chargino and neutralino has also been presented by the Compact Muon Solenoid 
(CMS) Collaboration I4 1 at the LHC. Therefore, high precision theoretical predictions for the 
associated production of xfx2 are very important for the forthcoming experiments at the 
LHC. 

The next-to-leading order (NLO) SUSY QCD corrections to the process pp xfxt 



mSUGRA was first investigated in Refs. [a, M where infrared singularities were dealt with 

n 

using the dipole subtraction method[7|. Also, the NLO SUSY QCD and SUSY electroweak 
(EW) corrections to this process in the general MSSM were calculated in Ref.|8|. In the 
following we further investigate the QCD effects on this process, including the NLO SUSY 
QCD corrections and, in addition, the next-to-leading-logarithmic (NLL) threshold resum- 
mation effects in mSUGRA using , hrtl at both the Tevatron 
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and the LHC 

The paper is organized as follows: In Sec. II we present the analytic results at fixed 
order. In Sec. Ill we briefiy summarize the threshold resummation formalism and derive the 
expressions for the resummed cross sections. In Sec. IV the numerical results are presented 
and discussed. Sec. V contains a brief summary of the conclusions. The SUSY vertexes 
involved in our calculations are summarized in Appendix |Al The abbreviations for the 
Passarino-Veltaman integrals are defined in Appendix [Bl The standard matrix elements 
and the explicit expressions for the form factors are summarized in Appendix O 

II. CALCULATIONS AT FIXED ORDER 

For hadron colliders the total cross section for the hadronic process, 

A + B^Xi+xl + X, (1) 

can be factorized into the convolution of the parton distribution functions and the parton 
cross section, 

Cr('S') = XI / fa/A{Xa, ^lf)fb/B{Xb, ^f)aab{s = XaXbS^Os), (2) 

a,b 

where /i/ is the factorization scale, f{x,^f) is the parton distribution function (PDF) and s 
is the parton center of mass energy. A and B both refer to protons at the LHC and proton 
and antiproton at the Tevatron, respectively. The parton cross section a is given by 

^ = YsjT.\^\"dPS'-\ (3) 

where ^ indicates the summation over final states and the average over initial states and 
J dPS^""^ represents the phase space integration. 

For simplicity, in this section we only present the expressions for the subprocess 

u + d~^xt + xi (4) 
The other processes are given by similar expressions. 

A. LEADING ORDER CALCULATION 

The LO Feynman diagrams are shown in FiglH Considering the light quarks as massless, 
at LO the production of X1X2 proceeds mainly via an s-channel exchange of a W boson, a 
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FIG. 1: The tree level Feynman diagrams for chargino and neutralino associated production. 



t-channel exchange of a down-type squark, and a u-channel exchange of an up-type squark. 
The LO amphtude is 

2 

Mo = Mso + J2(^to + Mto), (5) 



fc=i 

with 



I^TmT ' (6) 



and 

■^uO - 7, 772 ' 



ii-M? 



where Dl = gw^udl \f^-, ^l, ^ij, (^2^ ^k2i '^k2-> ^k2-> ^^e coefficients appearing 

in the SUSY couphngs and their exphcit expressions are given in Appendix El The standard 
matrix elements M'^ are given in Appendix O The LO amphtude and all of the NLO 
calculations in this paper are carried out in t'Hooft-Feynman gauge, s, t and u are the 
Mandelstam variables defined as 

s = (Pl+P2)^ i={.Pi-Pzf^ u={pi-pif. (9) 

In order to simplify the expressions, we further introduce the following modified Mandelstam 
variables: 

i' = i- Ml+, u' = u- Ml+. (10) 

After the D-dimensional phase space integration, the LO parton differential cross sections 
are given by 



dt'du' lQiTpT{l-e) \t'u' - sMl^ 1 ^ xt^ ^ ^ xi J ^^^^ 



5is + i + u- M|+ - Mlo)J2\Mo\^ 



where e = (4 — D)/2 and the function is the Heaviside step function. The exphcit 
expression for ^|A1oP is 

^ 1 [ ' 2sM~.M^oaia%Cl,Ct. ' [jai? + jbincmi - Mg,)(t - M|) 



6 (^-^i) 



4DLA^(t-M2+)(t-M2o) 2 J 4Di^L(M-M2+)(M-M2o) 2 . 



E ^siyu ^ x{ ' ' xi_ sr^ 



{i - Ml){i - Ml) {u - Ml){u - Ml) 

^AlDlit - Ml^){i - M|) + AAlDI{u - Ml^){u - M|) + SAlArDIsM~-,M^o 
^ ' ' s-M^^^ 

(12) 

The LO total cross section at the hadron colliders is obtained by convoluting the parton 
cross section with the PDFs in the hadrons A and B: 



j dxidx2 [fu/A{xi, /U/)/j/B(a;2, /i/) + {A ^ B)]a^ , (13) 

where is the Born cross section for ud —>■ XiX2- Obviously, the LO results are finite and 
free of singularities. 



B. NEXT-TO-LEADING ORDER CALCULATION 

The NLO QCD (including SUSY QCD) corrections for the production of X1X2 consist 
of the virtual corrections, generated by loop diagrams of colored particles, and the real 
corrections with the radiation of a real gluon or a massless (anti) quark. For both virtual 
and real corrections, we will first give the results in the dimensional regularization scheme 
(DREG) lli, in which, to restore supersymmetry, we modify the Yukawa coupling at the 
one loop level fl, S H Q : 

9w{qqx) = 9w{i - ^)- (14) 
ovr 

We will show the results in the dimensional reduction scheme (DRED) [ij] and compare the 
two schemes. 
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1. VIRTUAL CORRECTIONS 



The Feynman diagrams for the virtual corrections are shown in Fi£rJ2l In the calculations 
of v.tua. co„.ect,o,. we used *e ccnpute. p.o,.a. pad^ge Fo^CalcM - generate 
the one loop amplitudes and the self energies. The unrenormalized amplitudes for the virtual 
corrections are given by 

24 

-^y = E E UScDVn + fsUSYVn)M:,\ (15) 
n=l a,b=L,R 

where the explicit expressions for the standard matrix elements M^'' and the form factors 
fqcDVn ^^"^ fsusYVn given in Appendix O The ultraviolet (UV) divergence in the 
amplitude for the QCD corrections can be expressed as 



M 



V 



k=l dfc k=l 



UV Air e 



and the UV divergence in the amplitude for the SUSY QCD corrections is 



M 



V 



UV 4:71 ei ^"t-M^- V ^ dr 9 

k=l dk r=l 

+ E -^'o^^V [ - E ^t^K- + AMI - 2u 

k=l "fc r=l 
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(17) 



where Cp = 4/3 and 5?- = -R^i-Rji — -R^2-^j2- -^'^ the 2x2 matrix shown below, and is 
defined to transform the squark q current eigenstates to the mass eigenstates: 

, Ql \ ~ cos 6a sinOg , , ^ 

: ' R'=\ h (18) 

qji I \ — sin 9q cos % 

with < % < vr, by convention. Correspondingly, the mass eigenstates M^^ and M^^ (with 
< Mq,j) are given by 

M~ 

I = RmliR^y, (19) 

M| 







+ 





+ 



V J 



+ 



g 



+ 



+ 



FIG. 2: The Feynman diagrams for the virtual corrections to chargino and neutrahno associated 
production. 



with 



M5 = I 



Ml = M| + Ml + M| cos 2/3 (J|^ - sin^ 
M| = M| + Ml + M| cos 2/3e, sin^ Ow, 



Ag — /itan/3. 



(20) 

(21) 
(22) 
(23) 
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Here M~ is the squark mass matrix. Mq ^ and Ag are soft SUSY breaking parameters and 
fi is the Higgsino mass parameter. and are the third component of the weak isospin 
and the electric charge of the quark q, respectively. 

In order to remove the UV divergences above we renormalized the wave functions of the 
(s) quarks and the masses of squarks, adopting the on-shell renormalization scheme jl^. And 
the squark mixing matrix must also be renormalized. Denoting M^^q, g^o and go the bare 
squark mass, the bare squark wave function, and the bare quark wave function, respectively, 
the relevant renormalization constants are then defined as 

Ml, = Ml + SMl (24) 

g,o = (1 + ^SZDq, + ^6Zlqr, (25) 



and 



go = (1 + IsZIPl + l5Z^nPR)Q- (26) 



After calculating the self energy diagrams in Figj2l we obtain the explicit expressions for 
the above renormalization constants: 



a ' 



5M|=-^{Ao(M|)+4M|(5o + i?i)(M|,0,M|)+4Ao(M|)+4M|i?i(M|,0,M: 
-E'5jAVo(M|)}, 

k=l 

(27) 

SZi = ^{[Bo + B^ + MliB', + i?l)](M|, 0, M|) + {B^ + M|i?;)(M|, 0, M|)}, (28) 
= 2n{Mf-Ml ) E-^?^'^t^o(M|), (.,r = l,2, s^r), (29) 

SZI = i^^^lriBo + 5i)(0, Ml Ml), (30) 

k=l 



and 



5Z'^ = E(42)'(5o + fii)(0, M|, M|), (31) 

k=l 



where B'- = dBi/dp^ and Aq and Bi are the one-point and two-point integrals 17(] , respec- 
tively. Since we will factorize the coUinear singularities into the parton densities, as will be 
discussed below, the MS scheme for the renormalization of the initial quark wave functions 



8 



should be used here. However, the initial quark renormalization constants have no finite 
terms, except in the SUSY QCD corrections which are irrelevant for the PDF's. There- 
fore, the on-shell renormalization scheme is equivalent to the MS scheme for initial quark 
renormalization. 

As for the renormalization of the squark mixing matrix, the counterterm for the squark 
mixing matrix is defined as 

R> ^R' + 6R^, (32) 

where the counterterm SR^ can be fixed by requiring that the counterterm SR cancels the 
antisymmetric part of the wave function corrections pj] . The squark mixing matrix R 
counterterm can be written as 

= \ Y.{^Z% - 5Zl)Rl. (33) 

k=l 

The corresponding counterterms for the virtual amplitudes are given by 
Mc =2 (^^L + ^Zt)Mso + ^SZl J2 ^to + 2 ^ _ 

s=l s=l dj, 

2 2 

+ i E + ^ri^)(^2"V2 + MiHi)c'^iszi + 5zi) 

s=l k=l ds Jfc 



2 ^ u-M? 

S = l S = l 



2 2 



Z E E^TtW + 7^M^^^^^^'<^ + M3''b"2)CUSZl + 5Z%) (34) 



s=l k=l "fc 



s=l k=l dj^ ds 

' ' 5Zl ^ 5Z% 



1 ^ ^ 

2A^A^ v\k2 3 k2 3 n^_^2 u-M? 

s=l k=l "fe 

and can be factorized: 

3 

•Mc = E E (fQCDCn + ftusYcJM:'. (35) 

n=l a,b=L,R 
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The explicit expressions for the form factors fqcDCn fsusYCn presented in Appendix 
[Cl In Eq. fl35|) . the UV divergences in QCD and SUSY QCD corrections are given by 



M 



c 



QCD 
UV 



asCp 1 
Att e 



M 



sO 



2 t + 2M? 2 



and 
Mc 



SUSY _asCFl 

UV 4:71 e 

2 



k=l 



(36) 



fc=i 



t - MJ 



X^5fXM|-4M|-£ + 3il// 



J2 StStMl -AMI-U + 3M? 



fc=l 
2 2 



(37) 



>^>^ Cj>(M3^^af, + M3^^6 | 



2 2 



A;2, 



s=l 



respectively. From Eqs.( |T6i) .(|T71). ( l36i) and ( 1371) . we obtain 



(A^y + A^c) 



k=l 



and 



(A1y + 7W 



_ a^C^ 3 

47r e 



J](M'o + -^«o) 



(38) 



(39) 



fc=i 



The UV divergences above cancel, as they must. The renormalized amplitude at one- loop 
order is UV convergent 

{Mv + Mc) =0, 
but it still contains infrared (IR) divergences: 

a, r(l - e' 



UV 



IR 47rr(l-2e) 
Mc 



^^^^^^ 



as r(l-e) 



IR 47rr(l-2e) 



and 



{Mv + Mc] 
Here = -Cp and AX = 



^ a. r(l-e) fAnf4Y A\ 



Mr 



(40) 

(41) 
(42) 

(43) 
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The 0{as) virtual corrections to the differential cross section can be expressed as 



di'du' 167rs2r(l - e) t'u' - sMl+ 

Xi 



(44) 



X 



5{s + i + u- M|+ - M|o)2i?e ^{Mc + Mv)M; 



Xl X2 ' 

The IR divergences in Eq. fl44l) include both the soft and coUinear divergences, which cancel 
after adding the real emission corrections and absorbing divergences into the redefinition of 



PDF'sllSl, as will be discussed below. 



2. REAL CORRECTIONS 



The Feynman diagrams for the real emission corrections are shown in Figsl3] and |H 




FIG. 3: The Feynman diagrams for the real corrections without squark resonances in chargino and 
neutralino associated production. 
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FIG. 4: The Feynman diagrams for the real corrections with squark resonances in chargino and 
neutrahno associated production. 



After calculating the relevant Feynman diagrams the amplitudes for the real gluon emis- 
sion process 

+ dip2) ^ xtiPs) + xUpa) + gip^) (45) 
and the real massless (anti) quark emission processes 

u{pi) + g{p2) xtiPs) + xliPi) + d{p5) (46) 



and 



can be written as 



and 



d{pi) + 9{P2) ^ xtiPz) + xliPi) + u{p^) (47) 



23 



Mrg = Y. E fRGuMh'n (48) 

n=l a,b=L,R 



21 



■Mm = J2 E fmnMi, (49) 

n=l a,b=L,R 

respectively. The explicit expressions for the form factors /j^^n ^^'^ Irqu ^^'^ standard 
matrix elements Mq^ and Mq^ in Eqs.f l48p and fHOl) are given in Appendix [Ci 

The phase space integration for the real corrections will produce soft and coUinear singu- 
larities, which can be conveniently isolated by slicing phase space into different regions using 
suitable cutoffs. We used the two-cutoff phase space slicing method j2o|, which introduces 
two arbitrarily small cutoffs, 5s and 5^ to decompose the three-body phase space into three 
regions. 



12 



The parton level cross section for real gluon emission o"^ contains both the soft and the 
coUinear singularities and, in general, can be written as 



a = a -\- a + cr , 



(50) 



where and a^^ are the contributions from the soft and the hard coUinear regions, re- 
spectively, and o"^*" is the hard noncollinear par The explicit forms are described below. 

In the soft limit the energy of the emitted gluon is small, i.e. < 6s\/§/2, and the 
squared amplitude X^I-^RGp simply be factorized into the squared Born amplitude 
times an eikonal factor ^eik'- 



(51) 



where the eikonal factor $ejfc is given by 



iPl ■P5){P2-P5)' 

The phase space in the soft limit also be factorizes: 
where dS is the integration over the phase space of the soft gluon which is given by [2 



(52) 



(53) 



dS 



\3-2e 



dE^Ef^'dn 



2-2e- 



2(2vr) ,0 

The parton level cross section in the soft region can then be expressed as 



(54) 



a' = {A'Kas^ll') J rfP5(2)^|^or J dS^,,k. (55) 
Using the approach in Ref.j2o|, after integration over the soft gluon phase space, Eq. flS^ 



becomes 



with 



a = a 



rri 



27rr(l - 2e) 



(56) 



Al = 2Cf, = -ACf In 5,, = ACp In^ 5,. (57) 

In the hard coUinear region, E^ > 6sV^/2 and —6cy/s < Ui2 = (pi,2 — PbY < 0, the 
emitted hard gluon is coUinear to one of the partons. As a consequence of the factorization 



theorems 



2l| , the squared amplitude for the gluon emission process ( H5|) can be factorized 



13 



into the product of the squared Born amphtude and the Altarelh-Parisi sphtting function 22| 

for u{d) u{d)g, 

\Mrg\ — > (47ra,/i^ )> A^o + • (58) 

^-^ \ ZUi ZU2 ) 

Here z denotes the fraction of the momentum carried by parton u[(£) u{d) with the emitted 
gluon carrying a fraction (1 — z) and Pij{z,e) are the unregulated sphtting functions in 
D = 4 — 2e dimensions for < 2 < 1, which are related to the usual Altarelli-Parisi splitting 
kernels [22] as follows: Pij{z,e) = Pij{z) + eP[j{z). Exphcitly 

Puu{z) = Pm{z) = Cf^^ + Cf^6{1 - z), (59) 

and 

PLi^) = P'd^) = -Cf{1 -z) + Cf^6{1 - z). (60) 

The three-body phase space can also be factorized in the coUinear limit and, for example, 
in the limit —6cS < Mi < it has the following form[2o|: 

dPS^'\ud^ xtxlg) dPS^\ud^ xtxl s' = ^^) ^Qjrt-,^ dzdu^i-(^ - 

(61) 

Here the two-body phase space is evaluated at a squared parton-parton energy zs. The 
three-body cross section in the hard coUinear region is then given by{2o| 



a, r(l - e) /47r/i2 



27rr(l-2e)V s 



( - ^)^c '[Puuiz,e)fu/A{xi/z)f^/B{x2) 



dz / I — z\ 

+ Pdd{z,e)fd/A{xi/z)fu/B{x2) + {A^ B)] — l J dxidx2, 



(62) 



where f{x) is a bare PDF. 

After subtracting the soft and collinear region of the phase space, the remaining hard non- 
coUinear part a^'^ is finite and can be numerically computed using Monte-Carlo integration 
techniques 23|]. The result and can be written in the form 

rfa^=^^|A^i.GprfM^'\ (63) 

where dPS is the hard noncollinear region of the three-body phase space. 

In addition to real gluon emission, other real emission corrections to the inclusive cross sec- 
tion for A + B xfx2 NLO involve the processes with an additional massless (anti)quark 
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in the final state. Since the contributions from real massless (anti) quark emission contain 
initial state collinear singularities, we also need to use the two-cutoff phase space slicing 
method [20|] to isolate these collinear divergences. But we only split the phase space into 



two regions since there are no soft divergences. Consequently, using the approach in Ref . [20| , 
the cross sections for the processes with an additional massless (anti) quark in the final state 
can be expressed as 



da 



add 



^ a^iga xtxl + X)[fg/Aixi)fa/Bix2) + (A ^ B)]dxidx2 

{a=u,d) 



+ da^^^" r(l -e) fATTfij 



27rr(l-2e)V s 



-)^c'[Pug{z,e)fg/A{xi/z) fd/B {X2 ) 



(64) 



dz / I — z\ 

+ fu/A{xi)Pdg{z,e)fg/B{x2/z) + {A^ ^^^~\.~^) ^^i^^2, 



where 



Pdg{z) 



2K + (1-^) 
-z(l-z). 



(65) 
(66) 



The first term in Eq. represents the noncoUinear cross section. The parton cross section 
0"*" can be written in the form 



a 



c 



2s 



(3) 



(67) 



7(3) 



where dPS is the three-body phase space in the noncoUinear region. The second term in 
Eq. (l64ll represents the collinear singular cross sections. 



3. MA SS FA CTORIZA TION AND NLO TO TA L CR OSS SECTIONS 



As mentioned above, after adding the renormalized virtual corrections and the real correc- 
tions the parton level cross sections still contain collinear divergences. These can be absorbed 
into a redefinition of the PDF's at NLO, using mass factorization 22]. In practice this means 
that first we convolute the parton cross sections with the bare PDF's fa/uix) {H = A,B) 
and then use the renormalized PDF's fa/nix, fif) to replace fa/nix). In the MS convention 
the scale-dependent PDF's fa/nix, fif) are given by [20 1 

'^dz 



fa/nix, flf) = fa/H{x) + Y^ (^-^^ 




27rrfl 



-P:^{z)fp/H{x/z). (68) 
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This replacement produces a collinear singular counterterm which, when combined with the 



hard collinear contributions,gives, as in Ref. 
collinear contribution: 



20|, the 0{as) expression for the remaining 



da 



+ 



2ttT{1 - 2e) V s 



fu/A{Xi,^If)fa/B{x2, fif) + fu/A{Xi, fif)fd/B{x2, /^/) 

B, 



where 



fu/A{xi,iJ,f)fd/B{x2, iJ,f) + {A ^ B)^ a dxidx2, 



= Ci.(41n5, + 3), A^'= = AfM— ), 



y 



with 



P.M = PaM in(5,^4) - K^iy)- 

y 



(69) 

(70) 
(71) 

(72) 



Finally, the NLO total cross section ioi A + B ^ X^X2 ^^e MS factorization scheme 



IS 



a 



NLO 



a,f3=u,d 



/ dxidx2 \ [fa/A{xi,i2f)fi3/B{x2, fif)] {a^ + a^ + a^ + a 



HC\ I |_ ^coll 



a=u,d 



+ 5Z h^^i^^2[fg/A{xi,fif)fa/B{x2,fif) + iA<-^B)]a^{ga~^xtxl^)- (73) 



Note that the expression above contains no sing ularities for 2A^ + = and 2AY + Af + 
Af = 0. 



4. ON-SHELL SUBTRACTION 

In the massless (anti) quark corrections there is resonance production of squarks, which 
actually corresponds to squark and gaugino production at the LO followed by squark decay 



to a gaugino and a quark, as shown in FigjH We used the method in Ref. 13| to subtract 
their contributions. For example, consider a representative process 

u + g ^u + xl, u-*xt + d, (74) 
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which is shown as the first Feynman diagram in FigjH Using the Breit-Wigner propagator 
— + imT), the squared resonance matrix elements can be expressed as 

' ' (Q'-Miy + Mirr ^^^^ 



where Q = {p~+ + p^) . After subtracting the contributions due to resonance production 
the squared matrix element is 

m im (76) 

This subtracted result avoids double counting and makes the numerical calculation more 
stable since the resonance peaks are subtracted before the phase space integration. The 
dependence on the squark widths will be discussed in Sec. IV. 



5. NLO TOTAL CROSS SECTIONS IN BOTH DREG AND DRED SCHEMES 

In our calculations we used the DREG scheme. However, this scheme is not appropriate 
for SUSY models because it violates supersymmetry. To restore supersymmetry we modified 
the Yukawa coupling at the one loop level as shown in Eq. ffT^ . 

The real corrections and NLO total cross sections in the DREG scheme have been given 
above. Next we show the corresponding results in the DRED scheme. The contributions 
from soft gluon emission remain the same, but in addition to the the modified Yukawa 
couplings those from hard coUinear gluon emission and massless (anti) quark emission are 
also different. These differences arise from the splitting functions and the PDF's. 

First, note the LO amplitude in the DREG scheme with modified Yukawa couplings 
(DREGM) is different from that in the DRED scheme: 

2 

47r 



k=l 



Here, and below, the LO amplitudes and cross sections in the right hand side of equations 
are all in 4-dimensions, and their Yukawa couplings are not modified. Calculating the virtual 
corrections in the DRED scheme one finds that 6Zl, dZ^, dZf- and (5M| remain the same 
as in the DREG scheme. Thus 

J^DREGM _ J^DRED ^ q_ (^g) 
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However, the unrenormalized amplitudes Aiy differ: 

^DREGM _ ^DRED ^ (79) 

From Eqs. fl77|) . flTHl) and fl79l) . one finds the following relations: 

(A^o + Mv + Mc)'''''''''' - {Mo + Mv + Mcf''^'' = -^-Mo, (80) 

Air 

(^B ^ ^V^DREGM _ ^^B ^ ^V^DRED ^ _^^B ^ ^^^2y /g^) 

Second, note the splitting functions in the DRED scheme have no dependence on e: 

P.,(z,6)^^^^ = P,,(z). (82) 
From Eqs. (l69l) and (!82l) . one finds 

^^coU^nuEGM _ ^^coU^OREn ^ ^ ^^^^^ ^/„M(x,/y, /.,)PL(l/)L7i.(^2, /i/) 



a 



+ 2^ / —fa/B{x2/y,fif)P^^iy)fu/A{xi,fif) 

+ {A^ B)^a^dxidx2 + C(a^). 



(83) 



Third, note the PDF's in the DRED and DREG schemes are related 

(^y 

/3 

Substituting into the formula for the Born cross section we obtain an additional difference 
at 0{as) arising from the PDF's: 

(B\DREGM (R\DRED (^s j ST^ f U, ,, \DRED r>i / ,w / \DRED 

+ E /'-/«/B(x2/y,/./)^^^^P^Jy)/„M(a:i,/.,)^^^^ 
+ {A^ B)^a^dxidx2. 

(85) 

Finally note that Eqs. (183|) and (185|) are very similar except for the limits on the integral 
over y. Substituting Eqs. (IHTj) . (I83l) and (|85|) into Eq. (I73|l . we obtain the following relations 

18 



for the NLO total cross sections in two schemes: 

+ —fa/B{x2/y,fif)P'da{y)fu/A{xi,fif) (§6) 

J zn 

Using the exphcit expressions, including the e dependece, for the splitting functions P', one 
finds 

Therefore, the NLO total cross sections in the two schemes are the same at NLO. 
III. THRESHOLD RESUMMATION 
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iere we briefly summarize the basic formalism for threshold resummation Refs. 25|, |26|, 
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3G, 



31 



32|. Using pair inclusive (PIM) kinematics the invariant mass differential 



cross section can be written as 



u{S,Q'^)='^j dz j dXadXbfa/A{Xa, l^f)fb/BiXb, IJ,f)5{z - —^^)uJabiz,Q 



Q 



a,b 

where 



X(iXi)S 



2 



da ^ d&ab Q ,„„s, 

and is the invariant mass of the chargino and neutralino. The differential cross section 
il^ab contains large logarithmic terms a^[ln^{l — z)/{l — z)]^, which come from the incomplete 
cancelation between real gluon emission and virtual gluon corrections. In the region z 1 
[Q"^ s) these large logarithms have to be resummed to all orders in a^. 

In order to calculate the hard-scattering function ujab we consider the IR regularized cross 
section for parton-parton scattering which factorizes: 

UJab = j dz j dXadXb(l)a/a{Xa, lJ'f)(pb/b{Xb, IJ'f)5{z - ^ ^ ^ )6jabiz , Q^) , (90) 

where (pa/a and (pb/b are the flavor diagonal parton distributions in partons. Using a Mellin 
transformation with respect to r the convolution in Eq. ( l90i) can be simplified as the product 

COab{N) = 4>a/a{N)4>b/b{N)~Uab{N), (91) 
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where 



^{N) = J 


f dr V;(r), 


(92) 


a/a{N) = j 


^ dx x^"Va/a(a;), 


(93) 


Vfe(^) = j 


^ dx x^'-^cpb/bix), 


(94) 



and 



UJ 



dz z^-^uj{z). 



(95) 



Here the large logarithmic terms in Cjab turn out to be a" In"* and z 1 corresponds to 
N oo. The next step is to resum the logarithms of A^. 

In order to separate the soft gluon effects from the short distance hard scattering we can 
factorize the differential cross section into the form 

Q 

ij 



{N) = ^,/,{N)^b/mY.^uS\ 



(96) 



where I, J are color indices, Hjj describes the short distance hard scattering and Sji is 
a soft gluon function associated with noncoUinear soft gluons. The explicit definitions of 
Hjj and Sji can be found in Ref. [25i]. The ip^s are the center-of-mass parton distribution 
functions in which the universal collinear singularities associated with the initial partons are 
absorbed. 

From Eq.(l9Tl) and Eq.( l96i) we have 

^a/a^b/b, 



^ab 



-Tr[HS]. 



(97) 



4>a/a4>b/b 

After resumming the terms with the N dependence we obtain the exponentiated differ- 
ential cross section in the space of moments 25|, |26|: 



~EXP 
^ab 



exp 



exp 



X exp 
xPexp 
xPexp 



J22d, 



as I ^/5(as(/^'^)) 

fir 



Q 

dfl' 



(98) 
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where da^ is a constant and its definition is given in Ref . [27|] . P and P denote path ordering 
in the same sense as the integration variable /i' and in the opposite sense, respectively. The 
first exponent in Eg. (198 p resums the collinear and soft gluon emission from initial partons 
in the hard scattering and is given in the modified minimal subtraction (MS) scheme by 



(99) 



with 



ma.) = C,[^ + ^^K{^)y (100) 

z/(/^) = 2C;(a,/7r)[l-ln(2i;(;^))], (102) 

^(/0 = (A-^)VN', (103) 
Pi=Pi^/2/S, (104) 

where f3i is the particle velocity, n is the axial gauge vector, is the number of colors, and 
nf is the flavor number of light quarks. Cj = Cp = {N^ — 1)/{2N() for initial quarks and 
Cf = Ca = Nc for initial gluons. In Eq. (!98l) 7^ is the anomalous dimension of ip and is given 
at one loop by 7^ = 3CirQ;s/47r for quarks and 7^ = PqC^s/t^ for gluons. The P function is 
defined as 

f^M=y-^=-f:M-r^'\ (105) 

^ n=0 ^ 



with 



/?o = (llCA-2n/)/12, (106) 
/3i = {17Cl - 5CAnf - 3CFnf)/2A. (107) 



in Eq. (!98!) is the soft anomalous dimension matrix 26l| . which can be derived from the 
eikonal diagrams as shown in FiglSl and is given by 



27r 



-[-2 In 2 - ln(f(„)t;(j)) + 2 - 2m]. (108) 



Next Eq.( l98|) at NLL can be written in the simplified form 

c5f,^^ = ^oCabM exp [X(iV, as)] , (109) 
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FIG. 5: The one-loop eikonal diagram for chargino and neutralino associated production. 



with 



ujo = J dr t'^-'luo, (110) 

X{N,as)=gi{X)\nN + g2{X), (111) 

^i(A) = ^[2X + {1- 2A) ln(l - 2A)], (112) 

g,iX) = ^[2A + ln(l - 2A) + ^ 111^(1 - 2A)] " ^[2A + ln(l - 2A)] 



+^[2A + ln(l - 2A)] In ^ - ^2A In %, (113) 



where cuq = dao/dQ'^ is the Born differential cross section, N = N bxpIje), Ie is the Euler 
constant, and A = /3o«s IniV/vr. 

The function Cab can be expanded as 

oo 

a.(«.)=i+$:(v)"^^5^^- ^11^) 

n=l 

By matching the moments of the NLO cross sectionjs^ we obtain the function Cab in 
Eq. fllOQp at the NLO. The contributions to the first term C^^^ = C'^pas/ir in the expansion 
above come from the constant terms in the moments of the differential cross section, which 
are primarily the coefficients of the 5(1 — z) terms in the differential cross sections. The 
other terms come from the Mellin transformations of the logarithms depending on z. 

As shown in Eqs.(l38l). ( l39|) and ( 140|) in Sec. II, the divergences from QCD and SUSY 
QCD corrections cancel each other. Therefore, combining the contributions from real gluon 
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corrections and PDF renormalization, the QCD contributions to C^^^ are given by 
- \M^^ + ^1^ + 7^^^^^"^^^ + ^^ 



41niV-21n^ 



+41n2Ar-41n^lnAr+ — } , (115) 



with 
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n=l a,b=L,R 

Here the terms of order 0{ln N/N) and (9(ln(Q^/yU^)/A^) are included. 

In order to more completely include the behavior of the full towers of logarithms 
Eq. (fT09|) is modified: 



28] 



cu^r^ = cI;oexp[C^;^(a,)]exp[X(iV,a,)]. (117) 



To obtain the physical cross section we perform the inverse Mellin transformation, 

-C+ioo 



uj{t) = — / dNr-^'CuiN). (118) 



where the "minimal prescription" is used 



32|. 



To improve the convergence of the integration in Eg. 01181) we adopt the methods in 
Ref. 2^- First, we rotate the contour by an angle (j) with respect to the real axis and 
parameterize it in the form 

AT = C + ^exp^^"^, (119) 

where the upper (lower) sign applies to the upper (lower) half plane < z < oo (oo > ;z > 0). 
Then we rewrite the inverse transformation convolution: 

^(t) = ^/ ^Arr-^El(^-l)/«(^)H(^-l)A(^)]7t^ (120) 

'^'^ •^'^N a,b ^ ' 

where Ctv represents the modified contour. The inverse Mellin transformation of uJab/{N — 

rf, 
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is well behaved near the region z ^ 1 due to the suppression by the factor 1/(A^ — 1)^. The 
inverse Mellin transformation of (A^ — l)fi{N) is then 



-i. / dNx-''{N-l)MN) = -^[^M^)]=H^)- (122) 
^TTi Jqj^ ax 

Finally, 

^(r) = V r - /' —J'a{x)M-)yabiz). (123) 

And after integrating over the invariant mass in the differential cross section and inserting 
the terms ignored in the Mellin transformation we obtain the resummed total cross section 

(124) 



„RES „NLO I 



^NLL ^NLL, 

a —a n.~ 01 



"==0 " \ da 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

In the numerical calculations the following SM input parameters were chosen 

Mt = 170.9GeV[10], a{Mz)-^ = 127.918, ^.(M^) = 0.1176, ^^^^^ 
Mw = 80.403GeV, Mz = 91.1876GeV. 



The masses of the lig^ 
at the two-loop level 



it quarks were neglected. The running QCD coupling was evaluated 



33[ | and the CTEQ6.5M PDF's |3J] were used to calculate the various 



cross sections, either at LO or at NLO. As for the renormalization and factorization scales, 



we chose /i^ = /i/ = Qy±yO = ^/ (py± + p^o)"^, unless specified otherwise. 

Using the program package SPhenojss!] the MSSM spectrum, including the widths of the 
squarks, was calculated in the mSUGRA scenario in which there are five input parameters: 
the ratio of Higgs-field vacuum expectation values (VEV's) tan f3, the common scalar mass 
mo, the common gaugino mass mi/2, the trilinear coupling Aq, and the sign of the Higgs 
mixing parameter fi. The value of Aq does not significantly affect our numerical results so 
we put Aq = and, based on the analysis in the literature [2], [^gI, focused on /i > 0. 
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TABLE I: The dependence of the chargino mass M~± anc 
quark mass Mt- The masses were calculated using SPheno 



the neutralino mass M~o on the top 

A2 



351 for tan/3 = 5, mo = 200GeV, and 



m 



1/2 



250GeV. 



Mt/GeV 
170.9 
176.1 
180.1 



M~±/GeV 
188.506 
190.609 
191.978 



M~o/GeV 
189.915 
191.779 
193.003 



Table [T] shows the dependence of the chargino mass M~± and the neutrahno mass M^o on 
top quark mass Mt. We see that the chargino and neutrahno masses depend shghtly on top 
quark mass due to the fact that they are calculated using the SUSY renormalization group 
evolution (RGE). The explicit expressions for the total cross sections for the associated 
production of a chargino and a neutralino are independent of Mt as shown in Sec. II and 
Sec. III. Thus the top quark mass Mt only enters in the SUSY RGE. 

TABLE II: The NLO total cross sections for different squark widths T{lj) using and not using 
on-shell subtraction, respectively. 

r(g)/ro(g) (r'^Lo{PP^xt^2)lpb 
2 0.734 

1 0.825 
0.5 1.006 



C^NLoiPP ^Xtx2)/Pb 

0.650 
0.653 
0.661 



Table [IT] shows the NLO total cross sections for X1X2 production at the LHC, using on- 
shell subtraction, {o-nlo{pP XiX^))^ or not, {o-'^LoiPP XiX2))j different squark 
widths r(g), assuming tan/3 = 5, tuq = 200GeV and mi/2 = 250GeV. The squark widths, 
which were calculated using SPheno 35l|. are ro(g). Table HTl shows that the variation in 
^NLoiPP ~^ xtx2) is about 26% while the variation in aNLoijPP X1X2) is only about 2%. 
Obviously, using on-shell subtraction reduces the dependence on the squark widths. 

To present the resummation effects we defined the following quantities: 

da NLO 
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These represent the threshold resummation effects relative to the NLO cross sections. 

We present the numerical results for both xtX2 ^^"^ X1X2 production at the LHC, but 
show only those for X1X2 production at the Tevatron since these cross sections are different 
at the LHC but the same at the Tevatron. 

In FigllSlwe chose xtX2 production at the LHC as an example to show that it is reasonable 
to use the two-cutoff phase space slicing method in the NLO calculations, i.e. the dependence 
of the NLO predictions on the arbitrary cutoffs 6s and 6c is indeed very weak, as shown in 
Ref . {20I . Here aother includes the contributions from the Born cross section and the virtual 
corrections, which are cutoff independent. Both the soft plus hard coUinear contributions 
and the hard noncollinear contributions depend strongly on the cutoffs. However, these two 
contributions in (o-^o/t + c^hardcoii + virtual and (Jhardnon-cou) nearly cancel, especially for the 
cutoff 6s between 5 x 10~^ and 10"'^, where the final results for a^LO are almost independent 
of the cutoffs and very near T.lpb. Therefore, we will take 6s = 10"'^ and 6c = 5^/100 in the 
numerical calculations below. 

Using the same parameters we reproduced the results in Ref. {sl as shown in FigJT6l which 
provides a check on our calculations. However, our results are not exactly the same as the 
results in Ref.|5j because the masses calculated using SPheno[35] are different from those in 
Ref. si- 

FigJTTl shows the total cross sections as a function of tan/?, assuming mi/2 = 150GeV, 
for mo = 200GeV and lOOOGeV. The general shapes of the cross sections are similar. The 
main difference is that the absolute values of the total cross sections are different. X1X2 
production at the LHC has the largest cross section. In general, the total cross sections at 
the LHC are a few pb while those at the Tevatron are hundreds of fb. Fig{T7|also shows that 
the total cross sections for large tan/5(> 10) are almost independent of tan/5 while those for 
small tan/3(< 10) decrease with the increasing tan/? especially for ttiq = 200GeV. We note 
that the contributions from the resummation effects do not change the shapes of the curves 
very much. 

With the same parameters as in FigJlTlthe resummation effects 6K are presented in FigJTSl 
as a function of tan j3 for mo = 200GeV and lOOOGeV. Note that 6K is almost independent 
of tan/? for large tan/?(> 10) and there are larger resummation effects for mo = lOOOGeV 
than for mg = 200GeV. However, 6K at the LHC decreases with the increasing tan jS for 
mo = 200GeV and tan/? < 10. FigJTSl shows that 6K at the LHC for X1X2 production is 
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larger than that for X1X2 production and SK at the Tevatron is larger than at the LHC For 
mo = lOOOGeV the resummation effects can reach about 4% at the LHC and about 4.7% at 
the Tevatron. 

FigHH shows the total cross sections as a function of mi/2 assuming rriQ = 200GeV and 
tan (3 = 5. As mi/2 varies from ISOGeV to 250GeV M^o increases from lOlGeV to 190GeV 
and Mu^ increases from 406GeV to 599GeV, respectively. And the total cross sections 
decrease rapidly with the increasing of mi/2- For example, when mi/2 > 240GeV the total 
cross sections are less than Ipb and lOOfb at the LHC and Tevatron, respectively. Note that 
the total cross section for xtx2 production at the LHC is the largest while X1X2 production 
at the Tevatron is the smallest. 

With the same parameters as in FigllHlthe resummation effects 6K are shown in FigJ^U] 
as a function of mi/2- At the LHC the resummation effects increase with the decreasing of 
mi/2, reaching 3.3% for mi/2 = ISOGeV. At the Tevatron 6K increases with the increasing 
mi/2, reaching 4.9% for mi/2 = 250GeV. We also find that the smallest resummation effects 
6K at the Tevatron for mi/2 = ISOGeV are about 3.8%, which is larger that at the LHC for 
all values of mi/2. 

Figj2T] shows the total cross sections as a function of rriQ assuming mi/2 = ISOGeV and 
tan/9 = S. As mo varies from lOOGeV to lOOOGeV, M~± increases from 96GeV to 116GeV 
and M^o increases from lOOGeV to 117GeV, respectively. The total cross sections decrease 
with the increasing mo for mo > 300GeV. However, note that the total cross section for 
X1X2 production at the LHC is independent of mo for mo < 300GeV. 

With the same parameters as in FigEH the resummation effects 6K are presented in 
Figl22] as a function of mo. The resummation effects increase at both the LHC and the 
Tevatron as mo increases. For mo = lOOOGeV the resummation effects reach 3.9% at the 
LHC and 4.7% at the Tevatron. The total cross sections increase rapidly with the increasing 
mo for mo < SOOGeV while they become independent of mo when mo > 900GeV. 

Figsj23] and EH show the total cross section for X1X2 production at the LHC and the 
Tevatron, respectively, as functions of the renormalization scale fir and the factorization 
scale fif, and for iir = l^f^ assuming mi/2 = 2S0GeV, mo = 200GeV and tan/5 = S. The 
dependence in the LO cross sections at both colliders is increased by the NLO corrections 
and the dependence is slightly decreased by the resummation effects. The /i/ dependence 
in the LO cross sections at the LHC (Tevatron) is decreased by the NLO corrections and is 
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further increased (decreased) by the resummation effects. However, setting Hf = Hr = fiscaie, 
the resummation effects reduce the scale dependence at NLO. In fact, from Figl23]it can be 
seen that the renormahzation/factorization scale dependence of the total cross sections at 
the LHC(Tevatron) is reduced to 5% (4%) with the threshold resummation from up to 7% 
(11%) at NLO. 

Figsj25] and [26] present the differential cross sections as a function of the invariant mass 
Q~±~o assuming mi/2 = ISOGeV and tan /9 = 5 for mo = 200GeV and lOOOGeV, respectively. 
We see that the maximum in the differential cross section occurs at about Q^i-o = 230GeV 
and 280GeV for mo = 200GeV and lOOOGeV, respectively, and the differential cross sections 
decrease rapidly with the increasing Qr;±~o. The NLO corrections change the shapes of 

Al A2 

the differential cross sections, especially for Q~±~o < 300GeV. The threshold resummation 

Al A2 

effects enhance the NLO differential cross sections more at moderate values of Qr±ro and 

Al X2 

much less so at low or high values of Q^±~o. 

Figj271 shows 6K(i as a function of the invariant mass. In general, after slightly decreasing, 
6Kd increases more rapidly for mo = 200GeV than that for mo = lOOOGeV. The resum- 
mation effects are significant for large invariant mass. For example, for mo = 200GeV, 
6Kd is larger than 18% and 35% at the LHC and Tevatron for (5r;±~o > SOOOGeV and 

Al A2 

(5~±~o > 1200GeV, respectively. However, in general, 6K is only a few percent as shown in 

Al X2 

Figsim M and M 

V. CONCLUSION 

In conclusion, we have calculated the QCD effects in the associated production of xfx2 
the MSSM within the mSUGRA scenario at both the Tevatron and the LHC, including the 
NLO SUSY QCD corrections and the NLL threshold resummation effects. Our results show 
that, compared to the NLO predictions, the threshold resummation effects can increase the 
total cross sections by 3.6% and 3.9% for the associated production of X1X2 X1X2 
the LHC, respectively, and 4.7% for the associated production of xfx2 Tevatron. In 

the invariant mass distributions the resummation effects are significant for large invariant 
mass. The renormahzation/factorization scale dependence of the total cross sections at the 
LHC (Tevatron) is reduced to 5% (4%) with threshold resummation from up to 7% (11%) 
at NLO. 
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APPENDIX A 



In this appendix we summarize 37| the SUSY vertexes involved in our calculations. 



1. The chargino-neutralino-W vertex is 



with 



— 9W[ ^ Vil^j2 



(A2) 



(A3) 



v/2 ' 

where gw = e/ sin6'iy, -Pl = (1 — 75)/2 and Pr = (1 + 75)/2. 9y[r is the weak mixing 
angle. Z is the neutralino mixing matrix while V and U are the chargino mixing 
matrixes. The chargino index is i{= 1,2) and the neutralino index is j{= 1,2,3,4). 
Also we define Al = A]^ and Ar = A^j^. 




-tYiMPL + A^Pn) 




-iYiMPi + A^Pn) 



FIG. 6: The Feynman rules for the chargino-neutralino-W vertex. 



2. The chargino-quark-squark vertex is 



(A4) 
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with 



Cu — gwyudUiiRgi, 



(A5) 
(A6) 



where Vud is the {u^d) component of the CKM matrix and W is the squark mixing 
matrix. s{= 1, 2) is the index of the relevant squarks in the mass eigenstates and 
i{= 1, 2) is the chargino index. We define = C^^ and Cy = Cy. 



u 




ds 
d 




u 



ds 
d 




-tCpPi 



Us 



FIG. 7: The Feynman rules for the chargino-quark-squark vertex. 



3. The neutrahno-quark-squark vertex is 



^x".f = -x'jia%PL + b%PR)qq: " qia%PR + blPL)x%, (A7) 



where 



and 



hi, = - V2gweqta.n6wRl2Zji 



(A8) 



(A9) 

Cq and III^ is the electric charge and the third component of the weak isospin of the 
left-handed quark q. 
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-i{b%PL + a%PR) 



{a%PL + h%PR) 



FIG. 8: The Feynman rules for the neutralino-quark-squark vertex. 



4. The squark-Higgs vertex is 

^5?'i^(G) = -D^UradlH- - DT<dr,H+ - D^^'uJlG- - D'S''u*JnG+, (AlO) 



where 



and 



D'T = 9wVudKiKiM'^f3)Mw/V2 



= -9wVudRliRiiCos{2(3)Mw/V2. 
m and n are the indices of the relevant squarks in the mass eigenstates. 



(All) 
(A12) 



Ur, 



H- 



G- 

\ 



Ur, 



H+ 

1^ 



dn 



G' 



FIG. 9: The Feynman rules for the squark-Higgs vertex. 



5. The squark-W vertex is 

C^'w = iD^^[dl{d^Um) - {d^dl)um]W^ + iD^^[u*^{d^dr,) - {d^u*jL]W+ , (Al3) 

31 



where 



\ P 



W: 



: -iD'{f {p + kY 



dn 



dn 



FIG. 10: The Feynman rules for the squark-W vertex, p and k are the four-momenta of Um and 
dn in direction of the charge flow, respectively. 



6. The chargino-neutralino-Higgs vertex is 

(A15) 



where 



(Cf)ii = fl'M/sin/5(f/jiZj3 - ^ (tan 61,4/ + Z^s)), 
(Cl )ii = sin/5(ViiZj4 + (tan 6*14/ ^ji + ^^2)), 



V^2 
7=2' 

V2' 



and 



We define = {CE)i2, Cf = {Ci()u, = ((7f and = (Cf )i2. 
7. The SUSY QCD sector of the four-squark vertex is 



where 



ca pa pa pa pa 
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(A16) 
(A17) 
(A18) 

(A19) 



(A20) 



(A21) 




FIG. 11: The Feynman rules for the chargino-neutrahno-Higgs vertex. 

a and /3 represent the flavors of the relevant squarks. Here j, k and / are the relevant 
squark indices, r, s, t and u are the color indices of the relevant squarks. 

\ u 



X 

s / \ 



■ —nrfi \T'^^ T"' -I- 7^" T"- <N« "s"^ A , 



FIG. 12: The Feynman rules for the SUSY QCD interaction of the four-squark vertex. 

8. The squark- gluon vertex is 

Here i and j are the indices of the relevant squarks in the mass eigenstates. r and s 
are the color indices of the relevant squarks. 

9. The quark-squark-gluino vertex is 

= -V2gsT:AURlPR - RlPL)r%s + T{RIPl - R%PR)qrq:,s]- (A23) 
Here i is the mass eigenstate index and s is the color index of the squark. 
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FIG. 13: The Feynman rules for the squark-gluon vertex, p and k are the relevant four-momenta 
of squark q in direction of the charge flow. 




s Qi 




s Qi 



9 9 

FIG. 14: The Feynman rules for the quark-squark-gluino vertex. 
APPENDIX B 

In this appendix, for simplicity, we introduce the following abbreviations for the 
Passarino-Veltman two-point integrals -Bj, three-point integrals Cj(j), and four-point inte- 



grals -Di(j), which are defined as in Ref. 171] except that we use internal masses squared as 
arguments: 





= B.{Ml,AMl), 


Bt 




Bt 


= bS,o,mI), 


Bf 


= S.(M|+,0,M2J, 


Bt 


= B,{u,0,Ml), 


Bl 


= 5,(M|,0,M|J, 


Bf 


= B,iO,MlMl), 


B^ 


= i?,(M2 ,0,Mj ) 


Bl 


= bS,o,mI), 


Bl 


= i?,(0,M|,MjJ, 


Bt 


= S,(m,M|,0), 
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5r = i?.(M|,0,M2^), 

Bf = B,iO,Ml,Ml), 
Bf = B,iO,Ml,Ml), 
5[ = 5,(MJ,0,M|), 
Bt = B,{Ml,0,Ml), 



Bl = B,iMlAMi: 
B^ = B,iMl,0,Ml 



a 



(0,5,0,0,0,0), 
(0,M2 ,t,0,0,Mj 



(0,M|,t,0,0,M£), 



(S,Mio,M^+,0,0,M4 
(S,M|,,M|,0,0,Mj 
(0,M|+,m,0,0,M2J, 
(0,M|o,M,0,0,M, 
(t,0,M2+,0,M|,M?^ 
(0,5,0,M|,M^^,MJ 



Ms/' 



2 n A/r2 n ;\/r2 

iu,Mlo,0,MlO,MiJ, 
(t,M|,0,M|,0,MjJ, 



:u,Mi+,0,M|,0,M^ 
;n,M|+,0,M|,0,MjJ, 



Ao)(0,5,M|+,t,0,M|o,0,0,0,M| 



XS' 



Ao)(0, S, M|, 0, M|^, 0, 0, 0, M-^J, 
Ao)(m, M|, s, 0, 0, M|+, M|, 0, Ml, Ml 



D 



D,^,)it, 0, s, M|, M|+, 0, 0, M|, M,^^, M| 



Many of the above functions contain so: 



Passarino-Veltman integrals can be reduced 38|] to the scalar functions Bq, Cq and Dq. 



t and/or coUinear singularities, but all the 



And the explicit expressions for these singular scalar functions have been calculated pre- 
viously in a dii 
LoopTools [isj]. 



viously in a different context 39|]. The remaining IR finite functions can be calculated by 
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APPENDIX C 



In this appendix we collect the explicit expressions for the nonzero form factors in 
Eqs.f lT^ - fl^ . The standard matrix elements in Eqs. ffT^ - fl5S]) for the subprocess 

U{pi) + d{p2) ^ xt{Pz) + X2(P4), (CI) 



are defined as follows: 



1\ /foh 




V2Pal''UlU4Pb1fiV3, 


Ti jab 
IVI2 




VlPaV3V2PbV4:, 


TiiTab 
M3 




V2PaV3U4PbUl, 






V2PaV4U3PbUi, 






V2PaV3U4Pb^2'^l^ 


n jab 




V2PaV3UAPb A i's'"! ) 


/I jab 
Mj 




V2PaV3U4Pb^3Ui, 


T\/Tab 




V2PaV4U3Pb^2'^lj 


l\/fab 

Mg 




U4PaUiV2Pb i>lV3, 






U3PaUlV2Pb i>lV4-, 


njab 




^2 -Pa Pll^V3U4Pb ^2ltJ.'^l 


71 /T n.b 




U4PaYUiV2PbPil^,V3, 


71 /Tab 




U3Pa'J^UiV2Pb^l1f_LV4, 


71 /TClb 




V2Pa i)ll^V4U3Pb i)2l^iUl 






V2Pa'y''V3U4PbAlfiUl, 


Afab 




V2Pal''V3U4PbltMUl, 






V2Pal^V4U3Pb'y^Ui, 






V2Pal''V4U3Pbi)21^lUl, 


Mtl 




V2Pal''YV3U4:PblfiluUl, 


M^t 




V2PaYYV4U3Pbl^lyUi, 


M^l 




V2PaUlU4PbV3, 






V2PaUiU4Pb^lVs, 






U4PaV3V2PbhUl, 






t;2 Pa A"l^4 A A^3, 



where a and b are the left-hand index L or right-hand index R, while Ui = u{pi) and 
Vi = v{pi) are the spinors of the particle with momentum pi. 
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The nonzero form factors in Eqs. (fT5|) - (l35|) are the following: 



u 



u 



/5^W2 = E gj^^f^^, {(t - Ml)Bl + 4tX^ + 4ti?J + (t - M|J [(5 + At + 

/5cw2 = E ^J^^^ {(^" - ^D^o + 4t5S + 4t5J + (t - Ml) [is + At + 
-M|+ - 2M|)C2^ - 2(S + ^ - M|,)(Co^ + + Ao(M|)} 

fQCDV. = E 3^(!^!^;;^ ((- - Mi)5o^ + 4^5o^ + AuB^ + (u - M^J (^^^ 

+2{u - M.^JCo^ + 2(^ - Mi,)C/; - (s + i)Cl + uC'o + Ml^Cf, - Mlf'^ 
+2{u - M^+)C( - 2sCf - 2iCf + 2M|+Cf + (3m - M^+)Ci + sCl + tCl 
+AuCl - Ml+Cl - 2M|C| - 2suD^Q + 2sM?^D^ - 2smDJ + 2sMlp\ 
^-2iuD\. - 2uMhD\. - 2iM? + 2MhM? + 2iuD\^ - 2uMhD\n 



? + 2M|oM| D?. + 2u2D^ + 2iuD\ - 2uMI+D\ - 2uMI,D\ 

Us 10 X2 Ms io Z Z z X2 ^ 



-2iMlpl - 2uMlDl + 2Ml-,Mlpl + 2MImId\ + 2MD^3 - 2mM|+D^3 
-2nM|oD^3 - 2tM?.^23 + 2Ml^^Mlpl, + 2M|m2^D^3 + 2^2/}^ _ 2snD^ 
+2tML)^ - 2uMl+Dl - 2uM^oDl + 2sMlpl - 2tM|^L>^ - 2uMlpl 
+2Ml^Mipl + 2Ml,Mlpl - 2{-i + + M|)(n - MDD^,) + ^(M^j) 

/5^W3 = E 3^|f ^ )2 ((- - + + ^uBt + iu- Ml) (4 



+2(« - Ml)C', + 2(^ - M|+)Co^ - (s + £)Co^ + uC', + Ml^C^ - M^fl 
+2{u - Ml+)C{ - 2sCl - 2iCf + 2M|+Cf + (3m - M^+)Ci + sCf + tCf 
+4MCf - M|+C| - 2M|oCf - 2smD^ + 2sMlpl - 2suD\ + 2sMlp\ 
+2iuD\^ - 2uMloD\^ - 2iMlp\^ + 2MloMlp\^ + 2tuDj3 - 2uM'1^,D\^ 
-2iMlp\^ + 2M^oMlD\.^ + 2m2D^ + 2tMi:)^ - 2uMI+DI - 2mM|o£'^ 
-2tM? ^5 - 2uM? L)^ + 2Ml+Ml D^. + 2M\Mj, D'i + 2MD5o - 2uMI+DI^ 
-2^M|D^3 - 2tM?.^23 + "iMl^Mlpl^ + 2M\,Mlp\^ + 2n^D^ - 25uD^ 
+2t7iL)^ - 2uMl^D\ - 2uMyp\ + 2sMlp\ - 2iMlp\ - 2uMlp\ 
^2Ml^Mlp\ + 2Ml,Mlpl - 2{-i + + M|)(n - M,2jZ}^3) + Ao(il/42^ 
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-2(£ - M? ) (^-C^ + sD;J + sD? + (M|+ - m)D^3 + (2s + t + - M|+ - Ml)!^^ 
+(s + i - Ml^)Dl^ + 2sD3« + (s + t)D3"3)) 

-2(t - M? ) (-Co^ + sD" + sDl + (M|+ - k)L>^3 + {2s + Uu- M|+ - M|o)L>« 
+(5 + i - Ml^)Dl^ + 2sDl + (s + £)L'3"3)) 

JQCDVb — 2^ + -^13 + -^23 + -^3 + "^33 J 

s=l 



fRR _ sr^ 

JQCDV6 — 



s=l 

fRL _ 2a^2^v-Pi2Q'3 
jQCDve - 2^ 



fRR _ \^ 
JQCDV7 — 2-^ 



_ ^ 2b%M^oC^yDl,a, 

JQCDV7 - 2^ ^ 

s=l 

fgcDvs — ~ E ^"^13 + -^3 + -^33) 

s=l 

_ A^2biM~^C(ja, 

fqCDYS — ~ (-^13 + -D3 + -D33) 

fQCDV9 = ~ E §^ ("^23 + -D3 + -D33) 

fqcDvw — E ^"^23 + -O3 + -^33) 

J^2hi2M^Cfjas 

fgCDVlQ = E ^ ^"^23 + -O3 + -^33) 

s=l 



JQCDVll — 2^ 



37r 

s=l 
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fRL _ Qs2^f-Pl2Q's 

Jqcdvu - ^ 
Jqcdvu - ^ 
Jqcdvu - 2^ ^ 

s=l 



fRR _ _ 
jQCDVn — / J 



fQCDV13 — 



2 ag,M-+C^D°3a, 
37r 

s=l 

2 biM~^C(jD^,,as 



s=l 



Stt 



37r 

s=l 

Jqcdvu - 2^ ^ 

s=l 



12"s 



jQCDV\h — 2^ 



2 6f2^~+C^f^J3«. 



JQCDV15 — 2^ 



s=l 

'2 af2M~+Ce.Dj3«, 



.=1 

Jqcdvw - 2^ ^ 
Jqcdvig - 2^ 

Jqcdvu - - 2^ 1^ 

s=l 

,LR ^ biM~^M^oCl,D^,,a, 



Jqcdvu — ^ 



Stt 

Jqcdvis — ~ 2-^ ^ 
Jqcdvis — 2-^ 



37r 

s=l 



fRR _ Y^ "s2'-^y-^oo"« 



Stt 

s=l 



^ ^" /^s p)b 

fRL _ Y^ "s2'-^y-^oo"s 

JQCDV19 — 2-^ ~ 



37r 
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^ nd r^s pta ^ 

fRL _ "s2'-^t/-^00"s 

JQCDV2Q - ^ 

JQCDV20 — ~ 



f, 



37r 

s=l 



m=l s=l ^ 



2 2 . „J 



m=l s=l ^ ^ ' 



fsusYvi — ^Q^/- *\ (2^-f'-^w*^oo + ctl2Cy-Doo(s — M-^)) 

fsusYvi = ~ X/ X/ Svrfs -^M^^ (2^l-D^J^*Coo + 0^2^*^-000 (s - M^)) 

m=l s=l ^ 
m=l s=l ^ (im'^ ^ rfs'^ 

[26^2 (t - Mj^) (^5^ + {s + 2i + u- M|+ - M|)Cr + (2s + 2t + 2m - 2M|+ 

2 

n=l 

+2M^M^(f - M|^)(Cr + C2'")i?ii?i)] } 

m=l s=l V Jm'^ ds' 

+CIj hai^ii - M? ) (bI + {s + 2i+u- Mh - M^o)C^ + (2s + 2i +2u- 2ML 

2 

-M|o)C2-) RiX, + hi, (4(<,< + Ri,Ri){iB{ + Ao(M|)) - E i^r"^Ao(M|J 

n=l 

+2M,M^o(t - M|^)(cr + cr)RLRi)] } 



&yv2 = E E '"wf!^^:f" (^o + tC^f + M|,C2^) 

m=l s=l ^ Jj'' 

/#^Vy2 = E E iS-Mtr ^^" + + ""k^^^ 

m=l s=l V ds' 

2 2 

5^ ^ 3n(u - M^\(u - M? ) i-^M,bl,M^o{u - M|J(C^ + C[)Rl,R%C^ 

m=l S=l ^ Wm/V Us J 

+2C2^^2^fiCf (^ - MiJ{B^ + {-§ -i-u + Ml + M'^^+ M|)C^ 

+{-s - i+M'^^ + M|)CO + b%{C^{i{M?B^ + uB'[){Rl,R% + Rl,R%) 
2 

- E Dr'^'AoiMl)) - 2M^M~^{u - M?J{C'o + Ct)C^Ri,R%)) 

n=l 



fRR 

JSUSYV2, 
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171=1 S=l '^m'^^ / 

2 

- J] Dr"^Ao(M?J) - 2M~,M~^^{u - M?J{C', + Cl)C^Ri,R^,)) 

n=l 



ISUSYV^ ^2^2^ 3^/^ - Mg ) ^ ° + ^9 + ) 

m=l s=l ^ "s' 



m=l s=l 

2Mr,aM%R: 



m=l s=l ^ H '\ W } 

+ {s - MI-){M^oAl - M~^^An){Cl + 2C^)D^^ + {s - M%.)a%M~^y{^ Dl{s - M^) 
+ {s - Ml.)ai^M^oC^yDt{s - M^)) 

m=l s=l ^ W) 

-{s - MI.){M~,Al - M^oAnWo + 2Q)^I^^ + al^M^+Cf (Df + D^,) 
x{s - Ml.){s - M^) - a%M^C^{Dl + 01 + Dl){s - Ml.){s - M^)) 

/I.Vv2i = E E ^(if^'^^i^^^-^"^^ - ^--)^^^ - - ^-))^o 

+{s -MiZ]{M^oAr. - M^^Anm + 2Cl)D^^ + - Ml.)a%M^^C]:i Dl{s - M^) 
+{s - Ml-)ai,M^Ct,Di{s - M^)) 



2 2 



-{§ - MI.){M~^Al - M^An)n + '^C'M^ + ai.M-^C^Df + Df) 

x{s- Ml-){s - M^) - a%M^oC^{Dl + + Dl){s - Ml.){s - M^)) 

m=l s=l ^ 

-2AR{Ci + C{ + Cl)D'^') 

/5"^5yv22 = E E ^^ffz^(<2^^^(^o + + + Dl){s - M^) 

m=l s=l ^ ^ ' 

-2AR{Ci + C\ + Ci)D^^) 



m=l s=l ^ ^ ' 

/I^5yy22 = E E sJg-MM' lCgv^;^!^ - Ml) - 2A,{C^, + CI + C^)D-) 
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2 O/T /?" 

m=l s=l 

^ ^ 2(1' 7?*^ 

m=l s=l 



m=l s=l 



^s2' 

37r 



m=l s=l 



-ai,M~^CUDt + Z^fi + 2Z}f2 + ^fa + ^2 + ^22 + ^23)) 

JSUSYV24 — 2-^ + -^11 + -^12 + 

m=l s=l 

fRL _ ^"'m2^v'^sJ^sl^ml ( rtd , rid , rid , rid \ 

ISUSYV24 — 2-^ + -^11 + -^12 + -^13 J 

m=l s=l 

/■M _ V V '^^ih^u^^^sA^i DC ^r)c ^r)c , d" ) 

JSUSYV24 — 2^ 2^ 2^ [J^l + -^11 + -^12 + -^13J 

m=l s=l 

^ ^ 9,^<^ ^, TDU T)d 

fLL _ \^ ^"m2'-^y"'s-"-s2-"-m2 ( rid , rid , rid _i_ nd \ 

JSUSYV24 — 2_^ 2^ ^ l-^l + -^11 + -^12 + -^13^ 

m=l s=l 

iQCDC2- L.^^^^_M2J2^^0 + ^l) 
2 4^2 C^Q, j,J 

Sqcdcs - - 2^ 37r(ii -M|)2^^o + 

g—l \ Us/ 

Sqcdcs - - 2^ 37r(^-M|)2^^o + ^1 j 
fRR ^ AuDLa 



ISUSYCl 



2 

fRL _ A/,DlQ:s -2 tdP I pd2 rjg , tdP pu2 , rj? prf2\ 

JsusYCi - 2^ 37r(s - ) ' -"-ife-^o + -Di-^ifc + ^i^ik) 



2(22 



f _ J \^ 0:50^2 / T>d T/d T/d T/d \ ( T/d nd nd nd \ A ( A/T'^ \ 

JSUSYC2 - Z^\2^ 2^ 3 it - M? Mi - M- '- ^1^1 ~ ^k2^m2)\^kl^sl - ^k2^s2)^0V^^dJ 

s=l \k=lm=l ^ cLi'^ ds' 



Cijasa'^2 



fe=i 



42 



+4Ao(M|) + Ao(M|J - J2{RiRi - Ri,RifA,{Ml 



2 N 

dk' 



2 2 



/; 



LL 

SUSYC2 



E EE 



k=l 



/ pd pd pd pd W pd pd pd pd \ 4 / ji ^' 



s=l I, k= 



2 



fc=l 



RR 

SUSYd 



37r(t - M| )2 

A;=l 



2 ^ 

dk' 



2 2 



E 



- -Rfc2-Rm2)(-Rfcl-Rri " -Rfc2-Rr2)^o(A^'^ 



2 ^ 
nfe; 



- Ml){B',R^l + BlRt + + RtBl) + 4Ao(M|) 



fc=i 



Y,(.RI,R:, - Rt,R%TAo{Ml) + AoiMl) + AMIB{ 



k=l 



fRL 

JSUSYC2, 



E 



2 2 



^ 1 37r(M -M| )(m-M| ) ^ 

0^2^ sCy 



XI '^-^^l-^™! ~ -Rfe2-Rm2)(-Rfcl-Rri - -Rfc2-R"2)^o(M|j 



1 m,=l 



.fc=l 



k=l 



The standard matrix elements in Eg. fl48l) for the subprocess 

U{pi) + d{p2) xHPs) + X2(P4) + ^(Ps), 

are defined as follows: 

^Gl = vTPaV3V^Pm[e^ ■ {Pl-P3)]{Tnmn, 

Mg'a = V^PaV3U,PbUT[e^ ■ {pi - p^W^Un, 

= U^PaV^V^^PtfuTiT^mn, 
M^l = U,PaUTv^Pb/''Vs{Tnmn, 

= vTPaVsVTPbfv^iTnmn, 

M^\ = v^P,v^U:,P,j^-u^{T- 

M^', = v:^PaV,u,p,4^uTiT^ 
Mh\ = v:^Pah<u,p,4^vs{T' 

M^\o = V^PaA<U,Pbfv3{T' 



(C2) 



Imni 

-IX \ 
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M^\, = V^PaYu^U^Ph^^Vsie^ ■ ipi+P2)]iT^)mn, 

M^\^ = U^PaUTv^PhfAMTlmn, 

M^\, = vTPaV3V','PbfAMTnmn, 

M^W = V^PaVsV^PbfhMTnmn, 

M^\, = U,PaUTv^Pa''hV3{Tnn.n, 

^G18 = V^PaV,U3Pb4''h<iTnmn, 

^G19 = v'2PaV,U3Pb4^MTiTnn.n, 

where x, m, and n are color indices for gluons, up quarks and down quarks, respectively. 
The nonzero form factors in Eq. fHHl) are the following: 
LR sr^ '^9sO'i2Cu 



fLR _ sr^ 

JRGl — 



JRGl — 



RL _ ST^ 2(y(saf2C*f 



fRL _ sr^ 

JrG2 — 2^ 



S + t + ti4 - M2 - M|o )(t24 - M? ) 



RG2 — 

s=l 

cRR _sr^ s2 I ^+ ^0 Us 



[s + t + t,^- M2 - M|o)(t24 - M? ) 



pRR _ \^ 

JRG3 — 2^ 



fRL _ sr^ 

JrG3 — 2^ 



J + t + tM - - M|)(ti4 - Ml) {Ml 
^ 2g,a%Ct. {sXi-Ml^^- Ml, + Ml ) 



u 



ys + t + tM - Ml^ - M^ojih, - Ml) {Ml - u) 
2D,gs{AnM^:- A,M,o) 
{M^-Ss,){s + i + ti4 - M2 - M|o ) 

Xl ^2 

,n 2D,gMLM~.-ARM^o;) 

RGA - I ^^2 
H 
2 



- su){s + t + tu- Ml^ - Ml,) 

Al A2 



fRR _ _ \^ 



gsK.C^yM- 



, V- f W + ^24 - - M|)(ti4 - M?^ 



2 



/iJG5 — 2-^ 
s=l 
2 

fLR _ _ "Sp 

JrG6 — / , 



gsa'i^C'yM- 



xt 



f M + t24 - - M|)(ti4 - 

gsajCljM^o 

{s + u + kA - M2 - M|o )(t - M? ) 

Xl X2 "s 
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.LL 9sbiC!,M^^ 

JRG6 2^ + + - Ml^ - Ml,){i - Mj ) 



s= 



rRR ^ tJs^s2^ul^I^+ 

JRG7 + i+i^^ _ ^2 _ ^2 )(^^^^ _ Mj ) 

JRG7 + £ + _ - M2 )(t24 - Mj ) 

Xl X2 "s 

.i^i. ^_Y^ 9sa%Ct.M^o 

JRGS 2^^(s + i + - Ml, - Ml,){u - Ml) 

fRL 9sb%kM^o 

jRGs + i + _ - M'^){u - Ml) 

rRR _ '^ArDlQs 
JrG9 



fRL 
JRG9 



fRR 
JRGIO 



fRL 
JRGIO 



fRR _ 
JRGU — 



fRL 
jRGU 



fRR 
JRG12 



fRL 
JRG12 



fRR _ 
JRG12, — 



fRL 
JrG13 





hi){s + i + iu - Ml, - 

Xl 


M|) 




su){s + i + iu- Ml,- 
2ARDLgs 


M|) 




■S34)(s + t + tl4-M|+ - 


-M|) 




-hA){s + i + ii4-M^+ - 
2ARDLgs 


-M|) 




-S34)(s + M + t24 

2ALDLg, 


--M|) 




- S34){s + U + i24- M^^ 

2ARDLg, 


--M|) 




- S34)(s + i + iu- M^+ - 
2ALDLgs 






-h4){s+i+iu-M^+ - 

2ARDLg, 


-M|) 




-S34){s+i+iu-M^+ 

2ALDLgs 





Xl '^2 



fRR _ _ gs^s2^V 

7HG14 A. + ^ + 4^ _ - Ml^^iiu - Ml) 

2 ^„ 



/•LR _ _ gsO's2'^V 

JRG14 A. (5 + ^ + 4^ _ ^2^ _ M|)(t,4 - M|J 



(S + M + t24 - M|+ - M|o)(t - M|) 
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JRGI, ^^{s + U + 44 - - M|)(t - Mj) 



LR _ 980^12^11 



fLR _ \^ 
JRG16 — 2^ 

s=l 
2 

/rgi6 — 



fRR 
JRG22 



fLR 
JRG22 



i-RR 
J RG23 



S=- 

2 



fRR _ 
JRG17 — 2^ 

s=l 

2 

fLR _ sr^ 



s=. 
2 



{s + U + t24 


-Ml+-M%)ii- 

9shi2Ctj 




{§ + U + £24 


-Ml+-M\)ii- 

9sK2Cv 




{s + u + £24 


-M|+-M|)(£i4 

9sO-^2^V 


-Ml) 


{s + U + ^24 


Xl X-2 

9sO-i2^U 




{s + i + iu - 


-M2 -M|o)(£24- 

Xl A2 

gACh 





5's«s2^Lf 



JRG18 — 

s=l 

rLL ^ 9sOl2^U 

JRG.S + £ + _ - M|) (44 - M|; 

2 

s=l 
2 

s=l 
2 

f-Ril _ 
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The standard matrix elements in Eq. fHQl) for the subprocesses 

u{pi) + g{p2) ^ xt{p^) + x^(P4) + rf(P5), (C3) 

and 

d{pi) + (7(^2) ^ xl(P3) + X^(P4) + m(P5), (C4) 

are defined as follows: 

= UsPaUTu^PbV,{Tnmn, 
Mi = U,PahV3U'',Pb4''uT{T^)mn, 
Mi = U.Pafv.u'^PbAuTiT^Un, 
Mi = U,P,Yv3U'',Pb-f,uT{rnmn, 
Mi = U^^P^V.U.P.f^.uTiT^mn, 
Mi = M^P,t;3M4A/" A<(T^)„^n, 
Mi = U.PaUTu'^PbfhMTnmn, 
Mi, = V^PaV,ulP,fi>^V,{T^)mn, 
Mi, = U.ParVsU'^Pbf hl,uT{T^)mn, 
Miu = V^PaV,U,PbV^{T^)mn, 
Mis = v['PaV3U4PbV^{Tnmn, 
Mi, = v",PahvTu4PbfVs{T^)^n, 
Mi, = vlPafv^U,Pbi)^V,{T^)^n, 
Mi, = y'lPaYvTuAPblMTlmn, 
Mi, = U,PaV^V^PbfhVs{T-)mn, 
Mi, = UsPaV^v'^PbfhMTnn.n, 
Mi, = v'^PaV.U.PbfhvTiTnmn, 
Mio = v'lPaV3U,Pb4^AvTiTnn.n, 
Mil = UAPa^VsV^Pb^^ hl^vTiT^mn, 

The nonzero form factors in Eq.( H9i) are the following: 

rRB '^QsK^C'v ( -[e^-{Pi-PA)]s-{e^-pi)iu + (e^-Pi)Ml [e- . (p, - p, - p,)] 

^^Qi 2-. _ Ml 1 s{i + U + SU- 2M2^ - M| + M?J s + u + - Ml^ - M| 
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^« ^ A 2gsa%C^ ( -[e^- (pi - ^4)]^ - (e^ ■ Pi)k, + (e^ ■ p^Ml _ [e^ ■ (pi - p3 - P4)] 
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FIG. 15: The dependence of the total cross sections for the associated production of X1X2 the 
LHC on the cutoff Sg, assuming mo = 200GeV, mi/2 = 150GeV, tan/3 = 5 and dc = (^s/lOO. 
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FIG. 16: The dependence of the total cross sections on the renormahzation/factorization scale with 
the same parameters chosen as in Fig. 2 of Ref.^j. 
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FIG. 17: The total cross sections as a function of tan/3 for the associated production of xfx2 
the two colhders assuming mi/2 = 150GeV, mo = 200GeV and lOOOGeV, = and ^> 0. 
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FIG. 18: SK, defined as SK = {a^^^ - a^^°)/a^^°, as a function of tan/? for the associated 
production of the two colhders assuming mi/2 = 150GeV, mo = 200GeV and lOOOGeV, 

^0 = and > 0. 
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FIG. 19: The total cross sections as a function of mi/2 for the associated production of X1X2 
the two colliders assuming mo = 200GeV, tan/3 = 5, Aq = and )U > 0. 
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FIG. 20: 5K, defined as 5K = {a^^^ - a^^°)/a^^°, as a function of mi/2 for the associated 
production of Xi5^ the two colhders assuming mo = 200GeV, tan/? = 5, = and > 0. 
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FIG. 21: The total cross sections as a function of mo for the associated production of xfx2 the 
two colhders assuming mi/2 = 150GeV, tan/? = 5, Aq = and n > 0. 
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FIG. 22: 5K, defined as 5K = {a^^^ - a^^^)/o^^^, as a function of mo for the associated 
production of x^X2 at the two colhders assuming mi/2 = 150GeV, tan/3 = Aq = and )U > 0. 
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FIG. 23: The dependence of the total cross sections for xtx2 production on the factorization 

scale(a), the renormaUzation scalc(b), and both scales equal(c) at the LHC assuming = 
200GeV, mo = 200GeV, tan/3 = 5, = and > 0. /xq = {m~± +m^o)/2. R = (t{h scale) /<^{l^o)- 
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FIG. 24: The dependence of the total cross sections for xtx2 production on the factorization 
scale(a), the renormalization scale(b), and both scales equal(c) at the Tevatron assuming mi/2 = 
200GeV, mo = 200GeV, tan/3 = 5, Aq = and /x > 0. Here /xq = (m~± + m^o)/2 and R = 

CriHscale)/<y{l^o)- 
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FIG. 25: The invariant mass differential cross sections for the associated production of xfx2 
the two colhders assuming mi/2 = 150GeV, tan/? = 5, mo = 200GeV, Aq = and /x > . 
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FIG. 26: The invariant mass differential cross sections for the associated production of xfx2 
the two colhders assuming mi/2 = 150GeV, tan/? = 5, mo = lOOOGeV, Aq = and fi > . 
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FIG. 27: The dependence oiSKa, defined as SKa = {da^^^ -da^^°)/da^^° , on the invariant 
for the associated production of X1X2 colhders assuming mi/2 = 150GeV, tan/3 

mo = 200GeV and lOOOGeV, Aq = and fi > . 
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